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Abstract 

The problem of joint downlink cell association (CA) and wireless backhaul bandwidth allocation (WBBA) in 
two-tier cellular heterogeneous networks (HetNets) is considered. Large-scale antenna array is implemented at the 
macro base station (BS), while the small cells within the macro cell range are single-antenna BSs and they rely on 
over-the-air links to the macro BS for backhauling. A sum logarithmic user rate maximization problem is investigated 
considering wireless backhauling constraints. A duplex and spectrum sharing scheme based on co-channel reverse 
time-division duplex (TDD) and dynamic soft frequency reuse (SFR) is proposed for interference management in 
two-tier HetNets with large-scale antenna arrays at the macro BS and wireless backhauling for small cells. Two 
in-band WBBA scenarios, namely, unified bandwidth allocation and per-small-cell bandwidth allocation scenarios, 
are investigated for joint CA-WBBA in the HetNet. A two-level hierarchical decomposition method for relaxed 
optimization is employed to solve the mixed-integer nonlinear program (MINLP). Solutions based on the General 
Algorithm Modeling System (GAMS) optimization solver and fast heuristics are also proposed for cell association in 
the per-small-cell WBBA scenario. It is shown that when all small cells have to use in-band wireless backhaul, the 
system load has more impact on both the sum log-rate and per-user rate performance than the number of small cells 
deployed within the macro cell range. The proposed joint CA-WBBA algorithms have an optimal load approximately 
equal to the size of the large-scale antenna array at the macro BS. The cell range expansion (CRE) strategy, which 
is an efficient cell association scheme for HetNets with perfect backhauling, is shown to be inefficient when in-band 
wireless backhauling for small cells comes into play. 
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I. Introduction 

In the upcoming 5G cellular standards, the radio access network (RAN) is going to experience a paradigm 
shift rather than an incremental upgrade of the 3G and LTE RAN ||T|. Conventional cell splitting for single¬ 
tier cellular architecture will not be suitable because of the extraordinarily complex radio network planning 
for interference management, and a huge demand for additional backhaul resources. The multi-tier het¬ 
erogeneous network (HetNet) architecture, incorporated with flexible backhaul connections and advanced 
physical layer technologies such as large-scale multi-user multiple-input-multiple-output (MIMO) antenna 
systems, is a promising solution [|2|. With a large number of small cells (e.g., pico and femto cells) de¬ 
ployed inside the macro cell range to meet the capacity and coverage demands, the network topology will 
see a substantial change. The newly incorporated technological and architectural enhancements will impose 
additional constraints to the system design. In particular, cell association (CA) and radio resource alloca¬ 
tion mechanisms will be fundamentally different, and they will have great impact on the performance of 5G 
cellular networks. 

Limitations in the backhauling systems for small cells has yet been adequately considered in previous 
studies of HetNets [0. With a very dense deployment of small cells, it is impractical to have wired con¬ 
nections at every small cell site. Instead, in-band wireless backhauling, which allows low-cost and fast 
deployment of a new cellular network, is a desirable solution. By implementing plug-and-play small cell 
equipments which can use in-band over-the-air links to macro base stations for backhauling, the operators 
may scale up their multi-tier 5G cellular networks in a short time. However, the use of in-band wireless 
backhauling gives rise to additional source of interference, which subsequently impacts the cell association 
and resource allocation strategies. It is therefore important to consider various backhauling scenarios in the 
context of designing cell association and resource allocation (or interference management) methods. For im¬ 
proved interference control under wireless backhauling, sophisticated technologies, e.g., large-scale MIMO, 
should be considered. 

This work studies the cell association problem in a two-tier HetNet where the small cell tier uses in-band 
wireless backhauling. Large-scale MIMO is considered at the macro base station (BS) to mitigate intra- 


and inter-cell interferences by exploiting the spatial degrees of freedom @|. Specifically, the macro base 
stations are equipped with large-scale antenna arrays and wired high-capacity backhaul connections. The 
single-antenna small cell tier, on the other hand, relies on in-band over-the-air large-scale MIMO links to 
the macro cell tier for backhauling. A macro BS therefore acts as a hub of the in-band wireless backhaul 
connections for small cells within its cell range. In such a system, we can either: i) dedicate certain amount 
of bandwidth for in-band wireless backhauling for small cells throughout the network, or ii) dynamically 
adjust the bandwidth allocation for backhauling locally at each small cell. Based on these two in-band 
small cell wireless backhaul bandwidth allocation (WBBA) strategies, the corresponding joint CA-WBBA 
scenarios for the two-tier HetNet are considered and analyzed. The main contributions of the paper can be 
summarized as follows. 

• For in-band wireless backhauling for small cells, two bandwidth allocation scenarios, namely, unified 
WBBA (u-WBBA) and per-small-cell WBBA (p-WBBA), are investigated for joint CA-WBBA in a 
two-tier HetNet with large-scale antenna arrays at the macro BS. 

• A duplex and spectrum sharing framework based on co-channel reverse time-division duplex (TDD) 
and dynamic soft frequency reuse (SFR) is proposed for interference management in the HetNet to 
facilitate the implementation of large-scale MIMO at the macro BS and in-band wireless backhauling 
for small cells. 

• For the u-WBBA scenario, a relaxed convex optimization approach is employed to solve the mixed- 
integer nonlinear program (MINLP) for the joint CA-WBBA problem for maximization of sum of loga¬ 
rithmic user rates. A distributed algorithm based on a two-level hierarchical decomposition framework 
is proposed. 

• The joint CA-WBBA under the p-WBBA scenario is shown to be a non-linear and non-convex MINLP 
problem. The distributed algorithm proposed for the u-WBBA scenarios is extended to the p-WBBA 
scenario by considering separate and iterative cell association and WBBA design. A fast heuristic 
algorithm of low-complexity is proposed to provide sub-optimal solutions for joint CA-WBBA in the 
p-WBBA scenario for online implementations. The General Algorithm Modeling System (GAMS) 
optimization tool is also employed to solve the sum log-rate maximization problem. 

The remainder of this paper is organized as follows. In Section!]!] some closely related works are reviewed. 
The system model under investigation is then described in Section [TTT1 In Section [IV] we study joint CA- 


WBBA for a two-tier HetNet under the u-WBBA scenario, and a distributed algorithm is proposed based 
on relaxation of the cell association variables and hierarchical decomposition. Section |V] first extends the 
distributed algorithm to the p-WBBA case and proposes a fast heuristic algorithm for macro base station 
offloading. Numerical results are presented in Section [VH before the paper is concluded in Section [VD] 

II. Related Work 

A. Cell Association 

The received signal-to-interference-plus-noise ratio (SINR) is an indication of radio signal quality and is 
directly related to error and outage performance. It is therefore a natural criterion for cell association and 
has been employed in most concurrent wireless systems. However, in multi-tier HetNets, when the BS load 
(i.e., number of associated mobile terminals) is taken into consideration, the SINR-based cell association is 
in general not optimal. The macro base stations in a multi-tier HetNet typically have much higher transmit 
power than the small cells. Load imbalance will then occur with SINR-based cell association Q. That is, 
the macro cell is likely to be over-loaded, while the small cells will have much lower load. 

As pointed out in 0, in a multi-tier HetNet, optimal cell association is coupled with scheduling and ra¬ 
dio resource allocation, which gives rise to a combinatorial optimization problem. Large power variations 
between different tiers make the problem even more involved. An effective solution is cell range expansion 
(CRE) of small cells through a biasing mechanism Even though determination of optimal bias is a 
challenging problem, the potential of the CRE for load balancing and system throughput improvement has 
been recognized. The almost blank subframe (ABS) for adaptive resource partitioning is a means to elim¬ 
inate inter-cell interference when CRE is used 0. The unbiased (macro) cells intentionally leave specific 
sub-frames blank so that the offloaded users associated with small cells are scheduled in these subframes 
to avoid inter-tier interference. A relaxed optimization approach was employed in |8J to convert the com¬ 
binatorial optimization for cell association into a convex problem, which is mathematically more tractable. 
A distributed algorithm, which avoids heavy overhead, was proposed through dual decomposition. In 0, 
the quality-of-service (QoS) provisioning of multi-tier HetNets was taken into consideration for cell asso¬ 
ciation design. A unified framework for outage minimization and rate maximization was investigated, and 
distributed cell association algorithms were proposed. Stochastic geometry is another tool, which has re¬ 
cently been used for modeling and analysis of cell association in multi-tier HetNets IfTOTl IfTITl . It models 




the locations of both the BSs and users as random point processes and analyzes the system performance 
statistically. 

B. Wireless Backhauling 

In order to satisfy the stringent delay requirements in the future 5G systems lfl2ll . it is important to jointly 
consider the designs of the access network and the backhaul network. A few recent works in the literature 
have investigated how limitations of the backhauling system for small cells affect interference management 
in multi-tier HetNets. Inter-tier interference in uplink of femtocell-macrocell overlay HetNets was studied 
in Ifl3l . where unreliable backhaul with varying capacity at femtocells was taken into consideration. Macro 
BS offloading problem for the uplink of a HetNet, where small cells rely on capacity-limited heterogeneous 
backhaul systems, was considered in |f3| - Based on a noncooperative game formulation of the problem, a 
distributed offloading algorithm employing a reinforcement learning approach was proposed. Macro cell 
users can therefore optimize their performance by being offloaded to neighbouring small cells. However, as 
pointed out in Q, in spite of a rich body of literature on multi-tier HetNets, there have not been sufficient 
studies on the impact of small cell backhauling solutions on radio resource management of HetNets. Few 
studies can be found which consider small cell backhauling in the context of designing cell association 
schemes for multi-tier HetNets. 

In-band wireless backhauling, which uses macro base stations as hubs of the backhaul system and connects 
the small cells to the hubs via in-band over-the-air links, would allow fast and cost-efficient implementations 
of large-scale 5G networks. To the best of our knowledge, the problem of cell association in HetNets 
considering in-band wireless backhauling constraints and large-scale antenna arrays at the macro BS, has 
not been investigated in the literature. 

C. Large-scale MIMO 

Large-scale MIMO, also known as massive MIMO, is an emerging multi-user BS technology that signifi¬ 
cantly improves radio spectral efficiency in a system with only single-antenna user devices lfl4l . Large-scale 
MIMO is considered as a disruptive technology for the future 5G cellular standards lfL5ll . In 5G multi¬ 
tier HetNets, large-scale MIMO can be used to effectively eliminate inter-tier interference. For example, 
a hotspot HetNet communication scenario was considered in lfl6l . and a spatial blanking strategy of the 
macro cells equipped with large-scale MIMO was studied to mitigate inter-tier interference. By employing 



Fig. 1 

A TWO-TIER HETNET WITH A SINGLE MACRO CELL AND SEVERAL SMALL CELLS AS WELL AS MOBILE USERS INSIDE THE MACRO CELL 

RANGE. 

a large-scale BS antenna array, the RF beam can be focused into ever-smaller mobile terminal (MT) spots 
and thus making the interference more “localized” lfl7l . The co-channel interference can be effectively 
eliminated, and remarkable improvement in spectral efficiency can be achieved. Consequently, we are able 
to realize more efficient joint transmission coordinated multipoint (JT CoMP) and more robust interference 
suppression. 

When small cells in the network have to rely on wireless links to the macro BS for backhauling, large-scale 
MIMO is a promising solution which not only improves service to macro cell users, but also provides good 
(in-band) wireless backhaul support to the small cells without requiring additional spectrum and hardware for 
a separate backhaul system. The “localized” interference pattern of large-scale MIMO has to be considered 
in cell association and wireless backhaul resource allocation design of future cellular systems. 

III. System Model 

A. Large-scale MIMO Macro Base Station and Small Cells Using Wireless Backhaul 

Consider a HetNet as shown in Fig. [Qwith a single macro BS equipped with a large array of Nr antennas, 
Ns small cells deployed within the macro cell range, and Ajy randomly located mobile terminals. For 
convenience, we denote the set of MTs by U and the set of small cells by S. So = S U {0} is the set of all 
BSs, where the index 0 is introduced for the macro BS. As in Ifl8l . we consider a large-scale antenna array at 
the macro BS having a beamforming group size N g -C Nr- We further assume that the beamforming group 



size N g is greater than Ns. On the other hand, the number of MTs Ajy within the range of the macro BS is 
much greater than N g . Resource sharing must be used when a large number of MTs are associated with the 
large-scale MIMO BS. These assumptions render the wireless channels between the macro BS and all the 
other nodes orthogonal. When all the small cells use their wireless links to the macro BS for backhauling, the 
cell association scheme must also take into consideration the joint resource allocation for user information 
and backhaul transmissions. With orthogonality of the wireless channels to and from the large-scale antenna 
array at the macro BS, high spectral efficiency is achieved with spatial-division multiplexing (SDM) and 
spatial-division multiple access (SDMA). 

In order to facilitate the uplink-pilot-based channel estimation of small-scale fading for the large-scale 
MIMO structure at the macro BS, the two-tier HetNet in Fig. |T|must operate in a TDD mode such that chan¬ 
nel reciprocity is guaranteed [H4ll ||4|| . The small cells use the wireless links to the macro BS for backhauling. 
Each mobile terminal selects one cell, either the macro or a small cell, to associate with. The objective of 
downlink cell association is to maximize sum of the logarithmic rates of the MTs such that a good balance 
between throughput maximization and fairness is achieved Ifl9ll . 

In this work, two in-band wireless backhauling strategies, namely, unified WBBA (u-WBBA) and per- 
small-cell WBBA (p-WBBA), are investigated. The former introduces a unified bandwidth allocation factor 
j3 e [0,1], which is the fraction of bandwidth allocated for wireless backhauling for all the small cells 
within a macro cell range. The p-WBBA scheme, on the other hand, determines (3j 6 [0,1], the fraction 
of bandwidth allocated for wireless backhauling, for each small cell j in S. In Sections HID and [TV] we 
investigate joint CA-WBBA in the two-tier HetNet with these two in-band wireless backhauling strategies. 

B. Wireless Channel Model 

With centralized massive antenna array implemented at the macro BS, the large-scale fading is constant 
during one cell association period and across a wide spectrum over all antenna elements. Conversely, the 
small-scale fading fluctuates fast enough relative to the cell association period so that it is averaged out in the 
channel measurements. We normalize the small-scale fading power and consider the average channel power 
gain of the macro-BS-to-MT channels, denoted by Hk . The channel model therefore only reflects large-scale 
fading, which is primarily determined by the locations of the communicating parties and the geographical 
characteristics of the radio propagation paths. We consider the non-line-of-sight (NLOS) path-loss model 
for urban macro BSs suggested in 3GPP TR 36.814 V9.0.0 and log-normal shadowing ll2()l . 


TABLE I 


Parameter values in the radio propagation path-loss model 


Parameter 

Value 

Carrier frequency 

3.5 GHz 

System bandwidth 

5 MHz 

Macro BS height 

25 m 

Small cell BS height 

3 m 

Average building height 

20 m 

Road width 

8m 

Noise power spectral density 

-174 dBm/Hz 

dB-spread of the log-normal shadowing 


for macro BS <jbs 

6 dB 

for small cells a sc 

4 dB 


By using the model parameter values given in Table HI the large-scale channel gain ///, in dB is specified 
as 

H k = 27.3 + 3.91 x 10 log 10 d k + Z(a 2 BS ), (1) 

where d k is the link distance in meters from the macro BS to the A: 11 ' MT, Z(cr BS ) is the log-normal shadowing 
term in dB, which is a zero-mean Gaussian random variable with variance a BS . Similarly, we use Gj and 
Lj k to characterize the wireless channel between the macro BS and the j th small cell and that between 
the j th small cell and the k th MT, respectively. Gj and L jk also account for large-scale fading only. The 
downlink wireless backhaul path-loss Gj for small cell j in dB is 

Gj = 24.6 + 3.91 x 10 log 10 dj + Z(a 2 BS ), (2) 

where dj represents the link distance in meters from the macro BS to the j th small cell. A NLOS urban 
micro BS propagation model suggested in lf20ll is adopted for downlink in small cells, which gives 

L j)k = 36.8 + 3.67 x 10 log 10 d j)k + Z(a 2 sc ), (3) 

where dj tk is the link distance in meters from small cell j to the k ih MT. As in 0]) and ©, Z (<rf c ) is the 
log-normal shadowing term in dB, which is zero-mean Gaussian with variance a 2 sc . 
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Fig. 2 

An EXAMPLE OF THE REVERSE-TDD TIME SLOT CONFIGURATION FOR MACRO BS AND SMALL CELLS IN A TWO-TIER HETNET. 


C. Interference Management Considerations for Wireless Backhauling of Small Cells 

As discussed in Section III-A1 the TDD mode for duplexing is employed by the two-tier HetNet to fa¬ 
cilitate the massive MIMO feature of the macro BS. Several TDD schemes for HetNets with large-scale 
antenna array macro BSs were studied in flTQ. The co-channel reverse-TDD scheme, which uses reversed 
uplink/downlink (UL/DL) time slot configurations for small cells and macro BSs, has a preferred inter-cell 
interference pattern and thus is a suitable model for this work. A simple example of the time slot configura¬ 
tions for macro and small cells operating in the reverse-TDD mode is shown in Fig. [2] 

The reverse-TDD scheme is a natural choice because of the use of in-band wireless backhauling of small 
cells. For convenience, we call the small cell’s functionality for communicating with its associated MTs as 
the service functionality, and the functionality for in-band wireless backhaul communications as the back¬ 
haul functionality. From the wireless backhaul’s standpoint, the small cells act as mobile users associated 
with the macro BS. During a downlink time slot designated by the macro BS, the small cells must be in re¬ 
ceiving mode for the backhauling functionality. When a small cell is in receiving mode, for backhauling, it 
should at the same time receive transmissions from the MTs, on frequency channels orthogonal to those used 
for wireless backhauling. Otherwise, significant self-interference will occur at the small cell. In this case, 
we see that this time slot has to be an uplink time slot of the small cell as a serving base station. Similarly, 
when a small cell is in transmitting mode, it simultaneously transmits to its associated MTs and the wireless 
backhaul hub (the macro BS) on orthogonal frequency channels. The time slot is therefore a downlink time 
slot to its own associated MTs, while it is an uplink time slot for the macro BS. Therefore, in a small cell 
using wireless backhauling, we have to reverse the DL/UL time slot configurations designated by the macro 
BS. In this work, we consider that all the time slots are equally allocated for uplink and downlink. All the 
small cells use the same DL/UL time slot configuration, which is the reversed version of that designated by 


the macro BS. 
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Fig. 3 

Cell edge inter-cell interference between adjacent cells in a TDD system when conflicting UL/DL time slot 

CONFIGURATION OCCURS. 


In cellular systems operating under TDD mode, severe co-channel interference may occur at cell edges 
if different time slot configurations are employed by adjacent cells. This can be shown by the following 
example, which is similar to the cross time slot scenario in TD-SCDMA networks ll22l . With the DL/UL 
assignments as shown in Fig. 3 for adjacent cells A and B, it is a DL time slot in Cell A and a UL time 
slot in Cell B in the fourth time slot. If there are cell edge users MT-1 (associated with Cell A) and MT-2 
(associated with Cell B), which are close to each other, MT-2’s uplink transmission will introduce significant 
co-channel interference to MT-l’s downlink reception in time slot number four. In our problem, where the 
small cells use cellular links to the macro BS for backhauling, the worst-case (adjacent cell) co-channel inter¬ 
cell interference occurs between a small cell and its associated MTs. That is, the small cell’s uplink backhaul 
transmission (to the macro BS) will interfere with its downlink communication with the associated MTs if 
they are assigned the same frequency band. Similarly, uplink transmissions of a small cell’s associated 
MTs will interfere with the macrocell-to-smallcell backhaul link at the small cell receiver. Because the two 
interfering communication links share a common node, i.e., the small cell BS, universal frequency reuse 
is out of the question. Soft frequency reuse (SFR) |24| is a good compromise between spectral efficiency 
and interference control in such circumstances. Specifically, in this work, we consider that the bandwidth 
allocated for wireless backhauling is not reused for user transmissions, either globally (in the u-WBBA case) 
or locally (in the p-WBBA case). 

Assume that the total frequency band can be divided into orthogonal sub-frequency-bands. In a multi-user 
MIMO system operating in the “large-scale MIMO” regime, i.e., N T is much greater than the number of 
simultaneous associated terminals, SDM within a macro cell range is considered. A beamforming group 





















size N g is available for each sub-frequency-band. Therefore, when the number of associated terminals for 


the large-scale MIMO macro BS is large, joint frequency-division multiple access (FDMA) and SDMA is 
employed. On the other hand, in a small cell, bandwidth allocation must be performed between communica¬ 
tion with its associated MTs and the wireless backhaul to avoid severe co-channel interference between the 
two functionalities. 

IV. Joint Cell Association and Unified Bandwidth Allocation for Wireless 


Backhauling of Small Cells 


A. Problem Formulation for Joint Optimization 

We consider a “large-scale MIMO” regime at the macro BS for its downlink transmissions. The two-tier 
HetNet in Fig. Q] operates under reverse-TDD mode. A binary indicator variable E {0,1} is introduced 
for cell association status of mobile terminal k eU with the j th cell (j E <S 0 ), where 1 indicates association. 
It is straightforward that must satisfy 



(4) 




According to the proposed reverse-TDD-SFR framework for interference management, downlink trans¬ 
missions to the MTs associated with the macro BS are only interfered by small cell users, which is suffi¬ 
ciently small to be ignored. The downlink of a small cell, on the other hand, is subject to interference from 
other small cells’ downlink transmissions. Therefore, we define the downlink SINR in the macro cell as 
7o,fc = r "x ti k ? where Pq is the transmit power of the macro BS. Similarly, the SINR in the wireless back- 

n A 

haul downlink for small cell j is jj = -jpp- The received SINR of MT k associated with small cell j is 
7 j :k = ;V(|+V J f' x)L lk ’ w here Pj (Pi) denotes the transmit power of small cell j (7). 

When multiple MTs are associated with a single-antenna small cell j E S, equal resource allocation 
among these MTs is optimal []8]|. For the u-WBBA scenario, the fraction of bandwidth allocated for wireless 
backhauling is (3 E [0,1]. The remaining bandwidth is used for user data transmission. The throughput of 
the j th small cell is therefore given by 



( 5 ) 






where the superscript “U” indicates the unified WBBA scenario. Because of the wireless backhaul constraint, 
the throughput R\j given by © must not exceed the achievable rate of the wireless communication link 
between the macro BS and the j th small cell, which will be analyzed in the following. 

In the “large-scale MIMO” regime, the macro BS has an antenna array size Nr much greater than the 
number of simultaneously served nodes (including both MTs and the small cells). In this work, we use the 
approximate model in fl8l and assume Nr N N g and Nj; > N g . A large number of MTs associated with 
the macro BS is served by resource sharing, e.g., time sharing. When a mobile terminal k is associated with 
the macro BS, its achievable rate f? 0 Jc in the “large-scale MIMO” regime is determined by both the channel 
coefficient II/., and the number of MTs associated with the massive MIMO macro BS lfl8l . Given the unified 
small cell bandwidth allocation factor f3, the achievable rate of the k th MT associated with the macro cell is 

m 

n u (1 — 0)N g ( N T -N g + 1 ^ 

R 0,k = X 0,k^ -- log 2 1 + --7o,fc , (6) 

l^k & u x o,k \ Ng J 

where equal resource sharing among MTs associated with the macro BS is assumed. Similarly, the capacity 
of the wireless backhaul downlink for small cell j is given by 


cf = /?io g2 (i + 


Njn — Ng 1 

n7 


-7 j 


(7) 


The downlink wireless backhaul constraint requires Rj < CJ such that all the downlink traffic of the 
associated MTs of small cell j in the current cell association period can be accommodated by its wireless 
backhaul. 

Note that as illustrated in Section IIII-Cl in order to facilitate the wireless backhaul design, a reverse- 
TDD scheme for duplexing is introduced for the two-tier HetNet. As a consequence, consider a time slot 
configuration as in Fig. [2l the information a small cell transmits to its associated MTs in a DL time slot of 
the small cell is that it received through the wireless backhaul from the macro BS in the previous time slot, 
which is a DL time slot of the macro BS. Because the radio propagation path-loss and large-scale fading are 
constant over several time slots, the inequality R 1 - < Cj 1 reflects the overall requirement on the downlink 
wireless backhaul for the entire cell association period containing a number of time slots. 

For convenience, we denote 

r j,k = log 2 (1 + 7 j,k), V(j, k) eS X U, 

Nr — N„ + 1 


1 + 


N n 


7o ,k i V/c G U, 


r 0 ,k = N g log 2 








and 


Cj = log 2 + Nt ^ 9 + 1 lj S j yj e 5, 

which can be considered as base-line rate expressions of the downlink channels that do not change with Xj^s 
or (3. The optimization objective, which is selected as the network’s sum log-rate, as discussed in Section 
III-Al to achieve proportional fairness, is given by 

'0--P) r 3,k] V- su 


rU = 5Z J2 x i- k log 

jeS 0 keu 


j&S 0 


( 8 ) 


. ThkeU X 3,k 

We denote X = {xj^\ j G So,k G Id} and have the following optimization problem PI for joint 
downlink CA-WBBA with unified wireless backhaul bandwidth allocation factor, which is a mixed-integer 
nonlinear programming (MINLP) problem. 


PI : maximize R U (B,X) 

/3,X V 7 

subject to Xj y k G {0,1}, V(j, k) G <S 0 x U\ (9) 

y. xj,k = i, vfc g u; (io) 

jG5 o 

0 < /? < 1; (11) 
RV < cf, Vj G 5. (12) 


For the u-WBBA scenario, the parameters R 1 - , Cj 7 , and /? r in PI are specified as in ©-([X]). 

It can be shown that the objective function in ® is concave, and the continuous variable 8 and the binary 
variables Xj ^s are separable in ®. However, the constraint in (fl2l) is a nonlinear non-convex coupling 
constraint. By fixing the bandwidth allocation factor 8, the constraint in (fl2l) becomes convex constraint on 
Xj t k s. Therefore, we will relax the cell association variables x :i js and consider a decomposition approach to 
solve the joint CA-WBBA problem with one unified WBBA factor /3 for all the small cells. 


B. Relaxed Optimization and Hierarchical Decomposition 

We relax the binary cell association indicators Xj t k s in problem PI, which gives 

0 < x jtk < 1, V(j, ()g5 0 x U. (13) 

Fractional cell association indicators can be interpreted as partial association with different cells in a cell 
association period. A hierarchical decomposition approach ||23| . which contains an upper level primal de¬ 
composition and a lower level dual composition, is considered to solve the relaxed optimization problem. 








By first conducting an upper level primal decomposition for the problem with relaxed association vari¬ 
ables, we obtain an inner problem for cell association, and an outer problem for unified wireless backhaul 
bandwidth allocation. 

Given a j3 value for unified WBBA, the optimization algorithm begins with the inner problem Pl.l 
formulated as 


maximize 

x 






subject to 

Xj,k > 0, V(j, k) G So x U; 

(14) 


Xj t k = i, g u, 

(15) 


j&So 



Rj (X; (3) - Cf (X; ^) < 0, Vj e 5. 

(16) 


where Rff s, defined by F&' = V, k log „ r i ,k —, are concave functions of x 7 - fc s. Because all con- 
straints in Pl.l are linear, the inner problem Pl.l is a convex optimization problem. 

Once the optimal solution X* = {x* k ] (j,k) e S 0 x U\ is obtained for the convex inner problem Pl.l 
parameterized by /?, it can be used in the following outer problem PI.2 for the unified wireless backhaul 
bandwidth allocation: 


PI.2 : maximize 

0 

Af ( ,log(l-/?) + ^Bf(X*) 


j&So 


subject to 

0<J3<1; 

(17) 


R?(j3;X*)~.CV(j3;X*)<0,jeS, 

(18) 


where R^' (X*) is the function value of R 1 -' evaluated at X*. In order to obtain the solution to the master 
problem PI for joint CA-WBBA with u-WBBA, the two sub-problems Pl.l and PI.2 are solved iteratively 
until convergence. We call the iterations for solving Pl.l and PI.2 for the primal decomposition the outer 
iterations, which are different from the inner iterations of the lower level dual decomposition algorithm to 
be discussed in more detail in Section ITlI-CI 

Note that R l - ; 's are constants in the outer problem PI.2. Maximizing the objective function of PI.2 with 
respect to /? is equivalent to maximizing log(l — /?) only, which is a concave and monotonically decreasing 
function of /3. Problem PI.2 then reduces to a feasibility problem whose solution is the smallest feasible 
value of (3 given the constraints ( IT71 ) and (fl8l) . By invoking the wireless backhaul constraint in (fl8l) . the 





u-WBBA factor /3 must satisfy 


3 > max 


Ekeu x l^j,k 


Vj G 5 


(19) 


y>.u + <-j) 

The solution to the outer problem PI.2, i.e., the smallest feasible 3 which satisfies the wireless backhaul 
constraint, is therefore given by (fT9l) at equality. Since the complexity of evaluating (fl9l) is low, and the 
number of small cells within a macro cell range is limited by the massive MIMO beamforming group size 
N g , which is typically not a very large number, globally optimal 3 for the outer problem PI.2 can be easily 
identified. The complexity of solving the master problem PI with the primal decomposition method lies 
mainly in solving the convex inner problem Pl.l. 

In order to solve the convex problem Pl.l for cell association in a distributed manner such that the need 
for complex inter-tier coordination is avoided, we consider a distributed algorithm considered in [[8|, in 
which a lower level dual decomposition will be conducted for the inner problem Pl.l. Auxiliary variables 
Kj = x j,k, Vj G cS 0 are introduced so that the inner problem Pl.l is equivalently given as 


Pl.lr : maximize 

x,{K j: jeS 0 } 

subject to 


EE x j,k log (r jik ) - E K > l °e( K i 

j£So k£U jSSo 

Xj,k > o, V(j, k) G So xW; 

Xj t k = 1, V(gW; 

jes 0 

Kj = ^2 x jt k, Vj e S 0 ; 

k&U 

0 < Kj < Nu, Vj G Sq] 


( 20 ) 




ft/ 


( 21 ) 

( 22 ) 

(23) 


By jointly considering the objective function in (l20l) and the coupling constraints in (12T1) and (1231) . we obtain 
the Lagrangian in (1241) associated with problem Pl.lr, in which we have // = [/i 0 ,//i, • • •, 3n s ] T an d 
v = [0,vi, ■ ■ ■ :Vn s ] t - 


£ U (K, //, u) 

= EE x j,k log (r jt k) ~ E A /io«W) + E Bj(Kj E^ + E ^^-(i-«E x j,kXj,k 

jG«So k£U o o k&A o _ k^jlA 


= E E x i* i> og(r j)fe ) - Hj - z/j(l - /0)r j)fe ] + Kj [nj - log(Kj) + VjpCj ]. 

k£U jGSo j€So 


( 24 ) 





The corresponding Lagrangian dual function is therefore given as 


9(1*, ”) = Skv ) + 9k(v, v)i 


( 25 ) 


where we define 


g k (li, v) = sup ^ x i,k (log(rj, fc ) - nj - Uj( 1 - P)r j}k ) 

Xj,k,jeSo j eSo 


(26a) 


subject to Xjj,. >0, Vj 6iS 0 ; 


J2 x i’ k = L 


j&So 


9k{h, v) = sup V Kj (fij - log (Kj) + UjpCj) 

Kj,j£S 0 jeSo 

subject to Kj < N[/, Vj G S 0 . 


(26b) 


The following dual problem for cell association minimizes the Lagrangian dual function u) with respect 
to )jl and v. 



Dl.l : minimize 


(27) 


Having only linear constraints in the convex problem Pl.lr, the Slater’s condition is straightforward and 
thus strong duality holds ll25l . The primal problem Pl.lr can therefore be equivalently solved by solving 
the dual problem Dl.l. 

C. Distributed Joint CA-WBBA Method With Unified WBBA 

Based on the hierarchical decomposition approach investigated in Section IIV-B1 we next present our 
method for distributed joint CA-WBBA in a two-tier HetNet employing the u-WBBA strategy. 

As illustrated in Section lTV-Bl solving the inner problem P 1.1 for cell association is equivalent to solving 
the Lagrangian dual problem Dl.l. According to (126al) . given //^s and lA' fi obtained from the f th inner 
iteration, the sub-problem g k {g^\ z/T) implies that each user k 6 U simply chooses the BS that offers the 
highest revenue Q ( fi = log (r jtk ) — — i/ r> (1 — P)fj,k- This mechanism for updating ;r ? ./,s is expressed as 


follows 



( 28 ) 








where 


je<s 0 L J 

Note that according to (l28l) . the algorithm only assigns binary values to the cell association indicator vari 


jk ] = argmax log(r i)fc ) - nf - uf\l - f3)r jyk , \/k G U. 


(29) 


ables Xj t k s without introducing any form of relaxation. The proposed hierarchical decomposition approach 
actually gives exact solution to the MINLP problem for joint CA-WBBA problem with unified WBBA. 


On the other hand, KjS can be updated locally at each BS by invoking the first order condition for (126bl) 


as 


Kj t+i) = min jexp(//^ + vf* fa - 1), iV(/} , Vj G S 0 . 


(30) 


It is noted that the Lagrangian dual function u) given in (1251) is not differentiable on // or u. A sub¬ 
gradient approach is used here to update the Lagrangian multipliers in the algorithm. Specifically, with 
carefully chosen step sizes S^t) and 5 u (t), which are usually quantities that diminish with the iteration 
count t, the Lagrangian multipliers are updated as 



(31) 


and 


(32) 


where the operator [•]+ picks the maximum of zero and its argument. Convergence of the sub-gradient 
solution to the dual problem Dl.l is guaranteed with appropriate step sizes S fl (t) and S u (t) |[8} ll26l By 
alternatively updating the primal and dual variables until convergence, we obtain the solution to the inner 
problem Pl.l for cell association, denoted by X*. 

The outer iterations in the master problem PI are conducted for problems Pl.l and PI.2, which are 
obtained from the upper level primal decomposition. Specifically, the cell association solution obtained 
from the inner iterations for problem Pl.l is fed back to the outer problem PI.2. Based on our analyses in 
Section ITlI-B l the unified WBBA factor is updated in the outer iteration as 



( 33 ) 


which is the j3 value to be used for problem Pl.l in the next outer iteration. 













We summarize the distributed method for joint CA-WBBA with u-WBBA based on our hierarchical de¬ 
composition solution in Algorithm [[] The operations and responses of different network entities, i.e., the 
macro BS, small cells, and MTs, are specified in the design. 

V. Joint CA-WBBA with Per-Small Cell Wireless Backhaul Bandwidth Allocation 

According to (fl9l) . the unified WBBA factor f3 studied in Section [[V] is chosen as the maximum of the 
smallest feasible per-small-cell WBBA ratios. Because different channel and load conditions are experienced 
by the small cells, the minimum WBBA ratio for each individual small cell may vary a lot. If the small cells 
within the macro cell range can adaptively use different WBBA factors, better spectral efficiency can be 
achieved. We are therefore motivated to investigate the p-WBBA scheme for joint CA-WBBA design in a 
two-tier HetNet. 

A. Problem Formulation 

With p-WBBA, we no longer use a unified and dedicated frequency band for wireless backhauling through¬ 
out the macro cell range. Interference management in the HetNet therefore becomes more involved. Similar 
to the problem formulation for the unified WBBA case in Section [IVl we define the small cell throughput 
Ft?, macro cell user rate R k)k and small cell DL backhaul capacity C? in the following, where the superscript 
“P” indicates per-small-cell WBBA. 

We dynamically allocate a fraction of the frequency band in a small cell j for its wireless backhaul. A per- 
small-cell frequency allocation factor [3 3 e [0,1] is introduced to denote the fraction of bandwidth allocated 
for wireless backhaul in the j th small cell. The remaining bandwidth is used for the small cell’s service 
functionality. The throughput of the j th small cell in the p-WBBA scenario is 



(34) 


Given the small cell bandwidth allocation factor / 3j for every j 6 S, if the sub-frequency-band assignment in 
each small cell is done in a random fashion, the average number of small cells associated with the macro BS 
on each sub-frequency-band for backhauling is N b = Yljes Pi- Consequently, similar to the unified WBBA 
case, the achievable rate of the k th MT in the macro BS under p-WBBA is given by 



( 35 ) 





Algorithm 1 Distributed method for joint CA-WBBA with unified WBBA 

Input: A feasible initial value of the unified WBBA factor (3. Large-scale channel coefficients ///, s, GjS, 

and Ljj,s for all j E S and k E U. 

Initialization: Calculate rj )k s, r 0 ,fcS, and CjS, for all j E S and k Eli. 

Output: Optimal cell association X* and WBBA factor (3*. 

Execution: 

while Convergence not achieved (outer iteration s) do 
while Convergence not achieved (inner iteration t) do 

- MT strategy (inner iteration): 
for Vfc G U do 

Calculate X (i+1 ) according to (l29l) and (1281) . 

end for 

- BS strategy (inner iteration): 
for Vj E So do 

Calculate Kf +P, s according to (l30l) . 

Update and according to ([3Tb and (l32l) . 

end for 
end while 

Updated cell association X* ( ' ,s+ I j as X'* +1) obtained at convergence of inner iterations. 

- Small cell strategy (outer iteration): 
for V. j E S do 

Calculate locally smallest feasible f3 as 

0 _ T. k& u^k + ' >r i,k 

P - E keU *% + 1 Hrj,k+c j y 

Feed back locally smallest feasible /3 to the macro BS. 

end for 

- Macro BS strategy (outer iteration): 

Designate the maximum of the locally smallest feasible /3s as /3( s+1 ). 

Pass 3 {s+] ) to the inner iterations for the CA sub-problem. 

end while 

Return X*( s+1 ) and /3*( s+1) at convergence of outer iterations as the optimal values X* and /3*. 






where we still assume equal resource sharing among MTs associated with the macro BS. The DL wireless 
backhaul capacity for small cell j under p-WBBA is 


0 . , Nt — N„ + 1 

Cj = Pj log 2 ( 1 + --7 j 


(36) 


The system’s sum of the logarithmic MT rates in case of p-WBBA is given by 

(1 - )r jik 


R 1 = 


E x °’ k log ( 

k&A V ^ k 


N a 


keu x o,k 


r 0 ,k 


EE^iogV 

1.^7/ 


'or J,‘ 
' %j,k 


3 r+£*r=£*f- < 37 > 

j£S j£S 0 


jes keu ^/keu j , 1 

Denote b = {/3j; j G 5} and the joint downlink CA-WBBA problem under per-small-cell WBBA is 
formulated as the following problem P2, which is also a MINLP problem. 


P2 : maximize Rf ( b, X) 

jeS 0 

subject to Xj y k G {0,1}, Vj G «S 0 , k G U: (38) 

E x 3,k = 1, Vfc G W; (39) 

jeso 

o <Pj< 1, Vj G 5; (40) 

Pj< E x J' fc > e (41) 

k£U 

Rf < Cf, Vj G 5. (42) 


Constraint (|4TT) is used to assign Pj = 0 when there is no MTs associated with small cell j, i.e., = 

0. 

The objective function of P2, which is specified in (l37l) . is nonlinear and non-convex. Moreover, ([37l) 
is not separable in the primal variables b and X as in the objective function of PI. Also, the wireless 
backhaul constraint in (l42l) is a nonlinear and non-convex coupling constraint. As we can observe from 
the problem formulation P2 and equations (I34l) - (l36l) . with large-scale MIMO and p-WBBA for small cell 
backhauling, the small cell specific bandwidth allocation factor Pj not only determines backhaul and user 
rates of small cell j, but also dynamically impacts the user rates in the macro BS. A decomposition approach, 
which decouples the design variables with unified WBBA in Section IIV-B1 is not directly applicable. The 
constraint in (l42l) is a nonlinear coupling constraint, which can be made linear if we conduct cell association 
and wireless backhaul bandwidth allocation separately. 

In this section, in addition to a distributed algorithm similar to that proposed in Section [TV] we investi¬ 
gate centralized fast heuristic algorithms for joint downlink CA-WBBA problem under p-WBBA. Directly 








solving the MINLP problem by the BONMIN (Basic Open-source Nonlinear Mixed INteger programming) 
solver If27l on the GAMS platform ll28l is also considered to provide benchmarks for the numerical results 
in Section IVTl 

B. Separate Cell Association and Per-Small Cell Bandwidth Allocation for Small Cell Backhauling 

By examining problem P2, we observe that with per-small-cell WBBA, the wireless backhaul constraint 
(l42l) is an active constraint so the optimum must be achieved at equality of (l42l) for all the small cells. Given 
the cell association indicators {x^; (j. k ) e S 0 x Id}, the optimal bandwidth allocation factor ff can thus 
be determined locally at each small cell by equating (1341) with (1361) . With distributed cell association, Xj^s 
can be determined at each MT, as in Section ITV-Cl Centralized cell association, on the other hand, relies on 
a central scheduler, which can be located at the macro BS or at the radio network controller (RNC). 

A straightforward solution to the joint downlink CA-WBBA problem is to consider the two tasks sepa¬ 
rately and iteratively in two algorithm stages, which is similar to the hierarchical decomposition approach 
investigated in Section [W] Given the cell association variables, which can be obtained by any simple CA 
strategy, e.g., SINR-based cell association, the optimal bandwidth allocation factors d 3 conditioned on the 
cell association can be determined accordingly with the knowledge of the wireless channel states as 

ft = £*«'** tr ** (43) 

X^k&U X jM\ r j,k + Cj) 

The small cell WBBA factors /3jS obtained from (l43l) are fed back to the cell association sub-problem 
to update the cell association scheme. The utility function R p to be maximized is concave given that 6 3 s 
are constant. Therefore, similar to the dual decomposition solution for sub-problem Pl.lr, a distributed 
cell association algorithm based on Lagrangian dual decomposition can be employed here to solve the cell 
association sub-problem distributively in each iteration of the master problem P2. Specifically, a Lagrangian 
similar to (l24l) is formulated for the CA sub-problem with per-small-cell WBBA 

£ P (X, = \ d i + MGa) - /£ - - 3 i) r j-C ■ X! ^ - lo g ( K i) + VjPjCj ], (44) 

k£l4 j GiSo j o 

where we have defined the per-small-cell WBBA related terms d 0 = log ^1 — '' j and d 3 = log(l — 

fjj) for j G S. The subsequent dual problem formulation is the same as Dl.l, except that g k fu, v) and 
gx(p, v) are slightly modified to include d :j s and use per-small-cell WBBA factors ffs instead of a unified 3. 






Specifically, we substitute (log(rj^) — fj,j — z/j(l — fi)rj,k) in (I26al) by (dj + log(r'j,k)~dj~ z^(l — (3j) r j,k), and 

simply replace (3 in (126bl) by / 3j . A distributed algorithm based on Lagrangian dual decomposition similar to 
that in Section ITV-CI can be used in the p-WBBA scenario for the cell association sub-problem. 

C. Macro BS Offloading and Small Cell Load Balancing 

The use of per-small-cell WBBA factors inevitably increases the amount of message passing for the dis¬ 
tributed algorithm discussed in Section lV-Bl The corresponding system overhead may be an issue in practical 
communication scenarios. Especially separating the CA and WBBA stages for problem P2 is not based on 
decoupling of CA and WBBA variables as in the upper-level primal decomposition for PI. The distributed 
iterative algorithm in Section IV-Bl is sub-optimal. In this subsection, we consider designing fast heuristic al¬ 
gorithm based on the idea of BS offloading which helps to achieve satisfactory performance. Such a scheme 
should provide a good compromise between performance and overhead/complexity. 

We start with the SINR-based cell association. Other cell association techniques can also be employed at 
the initial cell association stage. In SINR-based cell association, each MT selects the cell which achieves 
the strongest SINR as determined in Section |IIIJ Then in each small cell, the wireless backhaul bandwidth 
allocation factor is calculated locally according to (1431) . With SINR-based cell association, the high power 
macro BS attracts more MTs and subsequently causes load imbalance between the two tiers. Some small 
cells may not make good use of their available wireless backhaul resources because of their small or even 
zero load, while some are overloaded due to heavy load or less favorable channel to the macro BS. We 
are therefore motivated to design a heuristic algorithm for joint CA-WBBA from the perspective of load 
balancing in a two-tier HetNet where large-scale MIMO at the macro BS and wireless backhauling for the 
small cells come into play. 

We first consider offloading MTs associated with the macro BS to small cells for sum log-rate maximiza¬ 
tion. A greedy algorithm is used to sequentially check each MT associated with the macro BS to see if an 
improved utility can be achieved by offloading it to a small cell. The details of the offloading algorithm are 
given in Algorithmic It can be shown that the complexity of this algorithm is of 0{NjjNs). 

By running the same offloading procedure for each non-empty small cell after the macro BS offloading, 
improvement in the utility R p is expected. An overall complexity of O(NuNg) will be introduced by this 
supplementary greedy algorithm. However, the numerical studies revealed that the improvement due to the 
extra small cell load balancing procedure is negligible. 












Algorithm 2 Offloading of the macro BS 

Input: SINR-based cell association X and the subsequent wireless backhaul resource allocation b. Large- 
scale channel SINR { 7 ^; j e S 0 ,k E U} and { 77 : j e 5}. 

Initialization: Calculate N b , R p , Rf,Vj E S and R p . 

Output: Offloaded cell association X*. 

Execution: 
for Vk eU do 

if ( x 0 ,fc = 1) then 
for j G 5 do 

Let X T = X and set £ 0 ,fc = 0, Xj^ = 1 in X r . 

Calculate new fij as ,3' using X 7 and (1431) . and N[ t — N h — (3j + /3'. 

Calculate new R p ' by updating R p ' and R p ' with 3', N' b , 70 ^, 7 j t k and 7 j. 
if Improved R p ' > R p is observed, then 

Update fa = B' p N b = N' h , R% = Fig, R p = Rf, R p = /?/" and X = X T . 
end if 
end for 
end if 
end for 

Return the updated X as X*. 

VI. Numerical Results and Discussion 

In this section, we examine the proposed joint downlink cell association and wireless backhaul bandwidth 
allocation schemes for the two-tier large-scale MIMO HetNet through numerical examples. In each simu¬ 
lation trial, we simulate one large-scale MIMO macro BS at the center of its cell range, which is a circle 
of radius 350 m. Ns small cells and Ny mobile terminals are randomly and uniformly placed in the cir¬ 
cular macro cell range. An example (random realization) of the network geometry for the simulation with 
Ns = 10 small cells and Ny = 100 mobile terminals is shown in Fig. HI 
The macro BS has a transmit power of 43 dBm (20 W) and the small cells have transmit power of 33 dBm 
(2W). In addition, we assume a 5 dB transmit and receive antenna gain for the small cells, as opposed to 
0 dB gain of the mobile terminals. The large-scale path-loss and shadowing model as described in Section 





III-Bl is employed. At the massive MIMO macro BS, an Nr = 100 antenna array with beamforming group 
size of N g = 20 is considered. Different Njj and Ns values are examined in the numerical examples for 
the proposed algorithms. The SINR-based cell association and the results for CRE of small cells flU are 
also evaluated under the small cell wireless backhaul constraint for comparison. The SINR bias used in the 
numerical examples for the CRE scheme is 3 dB. For each set of system settings, the numerical results are 
obtained by averaging over a sufficiently large number of simulation trials. 



Fig. 4 

An example of the network geometry used for the 

SIMULATION. 


Fig. 5 

Empirical CDFs of the MT rate for different cell 
association schemes under the u-WBBA scenario. 


A. Joint CA-WBBA With Unified, WBBA Factor 

The empirical cumulative distribution function (CDF) of the long-term MT rate in the two-tier HetNet 
employing the proposed distributed joint CA-WBBA algorithm for the u-WBBA scenario is shown in Fig. 
|5l The CDFs of downlink transmission rates for MTs with the SINR-based cell association and CRE schemes 
are also plotted. In order to study the impact of the small cell density and the density of MTs on distribution 
of the transmission rates for the MTs, the CDFs for different Ns and Nu values are considered. 

It is observed from Fig. [5] that under the u-WBBA scenario, the proposed distributed joint CA-WBBA 
algorithm has the CDF curves significantly shifted towards the right, compared with the SINR-based cell 
association, which is an indication of improved MT rate for most of the MTs in the network. To better 
























Fig. 6 


The Ro.a and R 0 .g MT rate values for various cell association schemes with different N s values and (a) N n = 100 

MTs AND (B) Nu = 200 MTs UNDER U-WBBA. 

demonstrate the improvement in overall MT rate achieved by the proposed joint CA-WBBA scheme, we 
define the MT rate at probability p, R p , in (1451) for p e [0,1] and show in Fig. [ 6 ] the values of R 0 .5 and R 0 9 
for different cell association strategies in different system settings. 

Pr{MT rate > R p } = p. (45) 

R 0 5 is the median MT rate, and R 0 9 is the rate that can be attained by 90% of the MTs in the network. Both 
of these parameters can be considered as indicators of the overall QoS provisioning of the two-tier HetNet. 

According to the QoS indicators i ? 0 .5 and f? 0 . 9 , the distributed joint CA-WBBA algorithm significantly 
outperforms the other cell association schemes under u-WBBA. Consider the SINR-based cell association 
as the baseline scheme. Up to 77% improvement in f ? 0 .5 and 87% improvement in i ? 0 .9 are achieved by the 
proposed algorithm with Ns = 10 and Nu = 100. The improvement in the QoS indicators by the proposed 
joint CA-WBBA algorithm is consistently around 60% in most of the scenarios examined. However, the 
CRE of small cells by SINR biasing provides negligible improvement over SINR-based cell association, 
which is obvious in both the CDF curves (in Fig. [5]) and the corresponding R p values (in Fig. [6]). Therefore, 
the CRE is not an efficient solution for cell association in presence of in-band wireless backhauling for the 
small cells. 

We then evaluate through simulations the performance of different cell association strategies in terms 
of the objective function value of PI. The sum of the logarithmic rates of all the MTs in the network is 
shown in Fig. [7] for different system settings. The superiority of the joint CA-WBBA algorithm based on 
hierarchical decomposition is still phenomenal, while applying the CRE strategy gives very small gain over 
the SINR-based cell association. The objective function value does not change much with Ns, i.e., the small 




















































Fig. 7 



MT rate (bits/channel use) 


Sum of the logarithmic of MT rates for the proposed 

DISTRIBUTED CELL ASSOCIATION ALGORITHM UNDER U-WBBA 
WITH DIFFERENT DENSITIES OF SMALL CELLS AND MTS. THE 
PERFORMANCE OF THE SINR-BASED AND CRE CELL 


Fig. 8 

Empirical CDFs of the MT rate in the two-tier HetNet 

FOR DIFFERENT CELL ASSOCIATION SCHEMES UNDER THE 

P-WBBA CONSTRAINT OF SMALL CELLS. 


ASSOCIATION SCHEMES ARE SHOWN FOR COMPARISON. 


cell density. This is a result of the dedicated in-band wireless backhaul design of the small cells, as all the 
small cells have to share the same resource pool for backhauling. Another interesting observation in Fig. 
[7] is that an optimal system load (i.e., number of MTs) in terms of the objective function value is located 
at about Nu = 100 MTs for the proposed joint CA-WBBA algorithm, regardless of the small cell density. 
On the other hand, the peak appears at around Nu = 50 for the SINR-based cell association and the CRE 
strategy. This implies that the proposed joint CA-WBBA algorithm is capable of taking full advantage of 
the multi-user diversity gain of the macro BS antenna array with Nt = 100 antenna elements, while the 
SINR-based and the CRE cell association schemes can only enjoy part of benefits of the spatial degrees of 
freedom offered by the macro BS antenna array. 

B. Joint CA-WBBA With Per-Small Cell WBBA Factors 

In this subsection we show numerical examples for the p-WBBA scenario. Performances of the proposed 
algorithms, including the distributed joint CA-WBBA algorithm based on decomposition methods and the 
centralized heuristic algorithm for macro BS offloading (OFL), are studied. The SINR-based cell association 
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Fig. 9 

The Ro. s and R 0 . 9 MT rate values of various cell association schemes with N s = 10 and N v = 100, under p-WBBA. 


and the CRE strategy are evaluated for comparison. The solution given by the GAMS BONMIN solver is 
provided as a benchmark. Performance of the small cell load balancing algorithm, which gives negligible 
improvement over the macro BS offloading algorithm, is not shown in the results. 

We show in Fig. [ 8 ] the empirical CDFs of the MT rate for different cell association schemes. In order 
to focus on the main findings and make the plot clear, we only show the results for Ns = 10 with two 
different Nu values (Nfj = 100 and Nu = 200). By exploiting the flexibility in spectral utilization of 
the per-small-cell WBBA strategy, the rate performance of the MTs is in general improved compared with 
the unified WBBA scenario. Among all the cell association schemes considered in the simulation for the 
p-WBBA scenario, the distributed joint CA-WBBA algorithm and the centralized heuristic algorithm for 
macro BS offloading achieved more desirable performance. Specifically, the joint CA-WBBA algorithm 
can better reduce the fraction of MTs in the very low rate region, while the macro BS offloading algorithm 
works best in increasing the probability that the MTs fall in the high rate region. This can be interpreted as 
a tradeoff between fairness and median MT rate. The CRE strategy again fails to achieve a significant gain 
over the SINR-based cell association. The GAMS BONMIN solver, although gives very similar CDF curves 
as the joint CA-WBBA algorithm in the medium to high rate region, sacrifices a relatively large fraction of 
the least favorable MTs by leaving them in the very low rate region. Therefore, it is less preferable in terms 
of fairness. Similar observations can be made from the comparison of the QoS indicators Rq _ 5 and R 0 .9 in 
Fig. 01 or Ns = 10 and Nu = 100. However, in the p-WBBA scenario, the difference between CA schemes 
in f ? 0 .5 and R 0 , 9 values is less significant compared with the results for the u-WBBA scenario in Fig. [ 6 l 

The sum log-rate achieved by the distributed joint CA-WBBA algorithm and the heuristic algorithm for 


macro BS offloading with different MT densities is shown in Figs. 1 0(a)] and |10(b)1 for Ns = 5 and Ns = 10, 
























































(a) (b) 

Fig. 10 

The Ro. 5 and R 0.9 MT rate values of various cell association schemes with Ns = 10 and Nu = 100, under p-WBBA. 


respectively, under p-WBBA. The average utility values achieved by the SINR-based cell association, the 
CRE strategy, and the GAMS BONMIN solution are shown for comparison. The GAMS BONMIN solution 
achieves the highest objective function value, which is the “best” from the optimization perspective. How¬ 
ever, because proportional fairness is considered in the optimization objective, the solution still has fairness 
issues, which coincides with observations from Figs. [8] and [9l The distributed joint CA-WBBA algorithm 
based on the decomposition method outperforms all the other schemes except the GAMS BONMIN solution 
in terms of the achieved objective function values. By comparing Fig. QHwith Fig. [ 7 ] it is worth noticing 
that employing the per-small-cell WBBA strategy improves the objective function value. Similar to the ob¬ 
servations from Fig. [71 here we also see peaks of the curves locating at around Nu = 100. The difference is 
that, under p-WBBA, the SINR-based cell association and the CRE strategy can also take full advantage of 
the multi-user diversity gain by having an optimal load approximately at Nu = Nt = 100. 

The main observations of Section [VI] are summarized in Table [TT] In general, the schemes under p-WBBA 
outperform those under u-WBBA, however with a higher system complexity. The GAMS BONMIN so¬ 
lution, which gives the best objective function value, fails to achieve good fairness among the MTs. The 
proposed cell association schemes for the p-WBBA scenario, on the other hand, achieve more balanced 
QoS. Specifically, the distributed joint CA-WBBA algorithm achieves the best fairness in terms of 90% MT 
rate, while the macro BS offloading algorithm gives the best median MT rate. 














TABLE II 


Summary of the observations on different CA schemes 


CA scheme 

Performance (R0.5, R0.9, sum log-rate) 

Complexity 

u-WBBA SINR CA 

Lowest in all u-WBBA schemes in all aspects 

very low 

CRE 

Slightly improved over SINR CA 

very low 

CA-WBBA 

Significant improvement over SINR CA and CRE; 

optimal with u-WBBA; in general poorer than p-WBBA 

medium 

p-WBBA SINR CA 

No significant difference from others in f?o .5 and R0.9', 

low sum log-rate 

very low 

CRE 

Not guaranteed improvement over SINR CA 

very low 

GAMS BONMIN 

Optimal in sum log-rate; good R0.5 but poorest R0.9 

highest 

CA-WBBA 

Best R0.9', good f?o .5 but no significant advantage 

over SINR CA; good sum log-rate, slightly lower than BONMIN 

medium high 

Macro BS OFL 

Best R 0 . 5 ; good f?o .9 but no significant advantage 

over SINR CA; good sum log-rate, slightly lower than CA-WBBA 

low 


VII. Conclusion 

We have studied joint downlink cell association and wireless backhaul bandwidth allocation in a two- 
tier HetNet where the macro BS is equipped with large-scale antenna array, and the small cells rely on 
in-band wireless links to the massive MIMO macro BS for backhauling. A reverse-TDD and soft frequency 
reuse framework has been introduced for interference management in the HetNet model. With the wireless 
backhaul constraint, the joint scheduling problem for maximization of sum of log rates for MTs has been 
shown to be a nonlinear mixed-integer programming problem. Two wireless backhaul bandwidth allocation 
scenarios, which consider global (unified WBBA) and local (per-small-cell WBBA) backhaul bandwidth 
allocation have been investigated. By showing a separable concave objective function in the unified WBBA 
scenario, a two-level hierarchical decomposition method for relaxed optimization has been considered, and 
a distributed joint CA-WBBA algorithm has been proposed. The algorithm has been extended to the per- 
small-cell WBBA scenario, which achieves improved spectral efficiency and network utility. Fast heuristic 
algorithms for cell association in the per-small-cell WBBA scenario has also been studied. According to the 
numerical results, performances of the cell association algorithms are not sensitive to the density of small 
cells in the two-tier HetNet. An optimal load (i.e., number of associated MTs), which is approximately equal 
to the number of antennas at the macro BS antenna array, has been observed. It has also been shown that the 
cell range expansion strategy fails to provide good performance when the wireless backhauling constraints 










for the small cells are considered for cell association. 
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